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Abstract 



The solar wind plasma from the Sun interacts with the 
Moon, generating a wake structure behind it, since the 
Moon is to a good approximation an insulator, has no 
intrinsic magnetic field and a very thin atmosphere. The 
lunar wake in simplified geometry has been simulated via 
a l^D electromagnetic particle-in-cell code, with high 
resolution in order to resolve the full phase space dynam- 
ics of both electrons and ions. The simulation begins 
immediately downstream of the moon, before the solar 
wind has infilled the wake region, then evolves in the 
solar wind rest frame. An ambipolar electric field and 
a potential well are generated by the electrons, which 
subsequently create a counter-streaming beam distribu- 
tion, causing a two-stream instability which confines the 
electrons. This also creates a number of electron phase 
space holes. Ion beams are accelerated into the wake 
by the ambipolar electric field, generating a two stream 
distribution with phase space mixing that is strongly in- 
fluenced by the potentials created by the electron two- 
stream instability. The simulations compare favourably 
with WIND observations. 



Introduction 

A wake structure is formed as the solar wind flows past 
the Moon. The Moon is to a good approximation an 
insulator, has no intrinsic magnetic field and a very thin 
atmosphere. Measurements by Lunar Prospector m and 
magnetohydrodynamic simulations pi of the magnetic 
fields on the Moon indicate a typical scale height of 
100km for miniature magnetospheres around the mag- 
netic anomalies (these crustal remnant fields will be ne- 
glected here). This results in the Moon being a sink for 
solar wind particles that collide with it. Directly be- 
hind the Moon, there will therefore be a void in which 
no solar wind particles are present. As the solar wind 
continues beyond the Moon, this wake is filled in. 

Explorer 35 and the Apollo sub-satellites made the 
first measurements of the lunar wake, albeit with very 
low resolution. They were able however, to detect a 
significant depletion in density behind the Moon. The 
WIND spacecraft used the Moon for a gravitational as- 
sist on December 27, 1994, passing at a distance of 6.5 



lunar radii (Rl), through the lunar wake. With all its 
instruments switched on, a slice of the structure of the 
lunar wake was obtained. These observations showed 
significant ion and electron depletions, rarefaction waves 
travelling away from the wake, increased electron tem- 
peratures, constant ion temperatures, counter-streaming 
ion beams, an increase in magnetic field strength and a 
number of electromagnetic and electrostatic waves [|l 0, 

Analytical descriptions of the processes involved in 
filling lunar wake have been considered 0, pi O . Ogilvie 
et al. p] detected the counter-streaming beam distribu- 
tion, creating a two-stream instability. Farrell et al. |q| 
then analytically considered a two-stream ion distribu- 
tion with a background of electrons as the process by 
which the wake is replenished. 

Simulations of this nature have previously been un- 
dertaken |1J, |l^ |lj]. Intriguingly, whereas Farrell et 
al. |12f| find that the ion dynamics dominate the inter- 
action, the higher phase space resolutions of Birch and 
Chapman |^, |lj| established that the electron dynam- 
ics dominate the evolution of the wake. In this paper we 
present new results from these simulations, including de- 
tailed particle density and electric field descriptions, a 
comprehensive analysis of the electron instabilities oc- 
curring in the wake, and a detailed description of the 
structure and dynamics of electron phase space holes 
generated in the wake that are not found in previous 
simulations. 



1 The Particle-in-Cell Simulation 

The simulation was performed using a self consistent, 
fully relativistic, coUisionless l^D particle-in-cell (PIC) 
code. In these simulations, the ion and electron distri- 
bution functions are represented by a collection of "su- 
perparticles" . These "superparticles" , along with the 
electric and magnetic fields, can exhibit both full ion 
and electron kinetics. The velocities of the "superpar- 
ticles" and the fields have vector components in three 
dimensions and all quantities depend on one spatial co- 
ordinate, Y, and time. Since V • B = 0, the component 
of B along the simulation is constant, but given this re- 
striction, the simulation is fully electromagnetic. The 
phase space resolution of the simulation scales as ^ y^. 
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Figure 1: Simulation geometry. 



Vconv = 25vthi and t is measured in plasma periods, w~^. 
This is equivalent to X = 0.2747?^^ where t is in plasma 
periods and R^ is the radius of the moon, i.e. half the 
width of the void. In reality, Vthe > Vconv > '^thi] this 
may affect the magnitude of the potential and effectively 
rescales X = Vconvt- 

The simulation run time is limited by the length of 
the simulation box. If a wake related structure passes 
over the boundary and interacts with other structures 
moving in the opposite direction, the simulation must 
be terminated since this introduces an unphysical distri- 
bution. In these simulations a 2048 grid cell simulation 
box has been used (grid cell = Debye length, \d) which 
equates to ±8i?L either side of a 256AC initial void. An 
initial void of 128AD has also been simulated to verify 
that structures found are not a result of the choice of 
the size of the initial void. 



where n is the number of particles per cell. As in Birch 
and Chapman ||l3[, ^ 2500 particles per cell were used, 
compared with ^ 80 particles per cell used by Farrell et 
al. ^. 

Figure |l| shows the geometry of the simulation. The 
simulation box begins directly behind the Moon, per- 
pendicular to the solar wind flow (—X direction) and 
parallel to the Interplanetary Magnetic Field (IMF) (Y 
direction). A higher dimensional simulation would be 
required to describe the lunar wake when the IMF is 
at oblique angles, and bulk flows other than in the so- 
lar wind flow direction, as the simulation is performed 
in the rest frame of the solar wind with zero convec- 
tion electric field. The IMF will pass almost unaltered 
through the Moon since the Moon can be assumed to be 
a perfect insulator and therefore no currents can flow to 
change the field. The simulation begins with a void in an 
otherwise Maxwellian distribution of ions and electrons 
which mimics the removal of solar wind particles by the 
Moon. Full shadowing is appropriate since the electron 
and ion gyroradii (1.2km and 47km respectively) are sig- 
nificantly smaller than one lunar radius ('--^ 1738km). As 
the simulation box moves with the solar wind, the simu- 
lation run time can be equated to a distance behind the 
moon, effectively building up a 2D picture of the wake, 
provided the IMF and solar wind flow remain constant. 

An ion to electron mass ratio of 20 has been cho- 
sen in order to evolve both ion and electron dynamics. 
This alters the simulation's ion plasma frequency and 
ion thermal velocity to be 
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since T^. — Ti and n^ — Ui. Since the actual solar wind 
velocity is 25 times the ion thermal velocity Vtuii the 
simulation also requires this relationship. Thus, the dis- 
tance behind the moon is given by, X = Vconvti where 



2 Results of the Simulations 

We begin by discussing the structures and instabilities 
found and comparing them with theoretical predictions 
and other simulations. The results will then be com- 
pared to the WIND observations at 6.5i?L- 

2.1 Particle Densities 

The density profiles along the simulation box [Y(i?i)] 
are plotted on a colour plot at regular intervals [along 
X(i?2^)], with colour representing density. Figure Q(a) 
shows the evolution of the electron density over the en- 
tire simulation. 

The first structure to evolve is the rarefaction wave 
travelling away from the wake at the ion sound speed, 
Cs- This is due to the electrons moving into the wake, 
and more electrons from the undisturbed plasma filling 
the space they leave. Ogilvie el al. H used analytic 
solutions of plasma expansion into a vacuum and showed 
that in such a step in density of a plasma, from A'q to a 
vacuum at x = and i = 0, the expression for particle 
densities at later times along x are given by the following 



Ns = No exp 
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where Ng is the density of particle species s, Nq is the 
ambient solar wind density and Cs is the ion sound 
speed: 
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where Tg is the electron temperature and rrii is the ion 
mass. These equations describe a rarefaction wave trav- 
elling away from the density step {x = 0) in the —x di- 
rection (away from the wake) at the ion sound speed and 
an exponential decrease in density in the +x direction 
(into the wake). In the electron density plot (Fig. ||(a)). 
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Figure 2: Distance along the simulation box Y{Rl) against time (equated to distance behind the moon) X{Rl) 
showing (a) electron density evolution (cm~'^); (b) ion density evolution (cm~'^); (c) charge density evolution 
(e/cm^); (d) electric field magnitude evolution (|E|/Vm^^). The simulation box extends to ±8i?i. The ambient 
solar wind electron and ion densities are Scm^"^. 



the rarefaction wave is approximately where the density 
has dropped to 4cm~'^. 

InitiaUy, the ions behave similarly to the electrons. 
Figure 0(b) shows the evolution of the ion density over 
the entire simulation. Again, there is a rarefaction wave 
in the ion density. This wave can be seen to be co- 
incident with the electron rarefaction wave, and there- 
fore also travelling at the ion sound speed in agreement 
with the above equations. By taking the electron density 
away from the ion density at all points in Figs, ^fa) and 
(b), a charge density plot can be produced (Fig. H(c)). 

Further behind the Moon, once a sufficient number 
of electrons have entered the wake, more structures be- 
gin to appear. After —X{Rl) « IORl enhancements of 
the electron density and ion density in the central wake 
region can be seen. These are due to the electrons and 
ions streaming into the wake and setting up a counter- 
streaming beam distribution. Such a distribution is un- 
stable and subsequently a two-stream instability grows. 
This instability is discussed in more detail in Sec. III-D. 
The obvious density enhancements emerging from this 
region are due to the electrons being trapped by electric 
potentials generated by the instability. Narrower density 
enhancements are also generated by the two-stream in- 
stability. In these structures, the electron density forms 
a double peaked structure about a single peak in the ion 
density. From Fig. 0(c), the charge density plot, it can 
be seen that this is the main deviation from quasineu- 
trality in the simulations. These are electron phase space 
holes coupled with propagating ion structures which we 
will see are ion acoustic solitons and are examined in 
more detail in Sec. III-E. 

The majority of the structures associated with the 
lunar wake are confined to within the rarefaction wave. 
Some electron phase space holes, however, are acceler- 
ated above the ion sound speed and overtake the rar- 
efaction wave. The other detectable structures again 
occur after —X{Rl) w IORl, generally just outside the 
wake. These are slight fluctuations in the electron den- 
sity, travelling away from the wake at velocities between 
~ l.Svthe and ~ 3Avthe- They are caused by the most 
energetic electrons which have travelled right across the 
wake and out on the other side, setting up a bump- 
on-tail distribution. This distribution is unstable and 
again, an instability is grown, resulting in these electron 
density fluctuations. This is an example of the elec- 
tron bump-on-tail instability and is discussed further in 
Sec. III-F. 

Comparing these density plots with those produced 
by the simulations with an initial void 128Ad wide shows 
very little difference. The main difference is the distance 
between the double electron density peaks. They appear 
wider than in the 256Xd simulations, however, they are 
actually similar sizes in Debyc lengths (no t shown here). 

The simulations by Farrell it et al. ||l2| , display some 
similarities with these higher resolution simulations. 
Rarefaction waves and a two-stream instability occur, 
which shall be discussed further in Sec. III-D. No elec- 
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Figure 3: Electric potentials (V), beginning at 6.90;^^ 



1.9i?L, then every 13. 8w 



3.8Rl. 



tron phase space holes are simulated nor is the bump- 
on-tail instability detected. 

2.2 Electric Fields and Potentials 

Figure ||(d) shows the magnitude of the electric field 
on the same X, Y scale as Fig. g(c). The first electric 
field structure to appear is the ambipolar electric field. 
Significant increases in the electric field magnitude can 
be seen at Y{Rl) = ±1Rl. The electric fields at these 
points are directed into the wake. Since the electrons 
have a higher thermal velocity than the ions (Te = Ti 
and TTT-e < rrii) they move into the wake first, setting up 
the electric fields pointing into the wake. These fields 
then accelerate the ions into the wake. These initial 
charge separations can also be seen at the beginning of 
the charge density plot. 

Figure ^ shows the electric potential calculated at reg- 
ular intervals from the electric field data. The ambipo- 
lar electric field generates the drop in electric potential. 
This eventually disappears as the ions are accelerated to 
cancel the charge separation developed by the electrons. 
The average electric field magnitude decreases back to 
almost zero, after which the two-stream instability be- 
gins to grow. 

The two-stream instability produces a number of 
peaks in the electric potential. These peaks are associ- 
ated with the bipolar electric field structures, which ap- 
pear as a double peaked structures in the electric field 
magnitude plot (Fig. 0(d)). Electrons are trapped in 
these electric potential peaks, oscillating back and forth. 
Some of these peaks evolve to become electron phase 
space holes and move out of the wake, coinciding with 
the double electron density peak about a single ion peak. 

The electric potential peak which grows wider as the 
simulation evolves, remains nearer the centre of the 
wake. This still has electrons oscillating back and forth 
as they are confined within the peak. In Fig. g(a), the 
electron density plot, oscillating density perturbations 




F 


-1 1 
V(R1 


"H 


m 




m 


w. 




-1 


-2 


2 -2 


-1 



MR) 


1 


2 






Figure 4: Electron phase space evolution [Y{Rl) against 
vy], beginning at Oujp^ => ORl, then every ll.ltjpg^ => 
S.ORl- 



Figure 5: Ion phase space evolution [Y{Rl) against 
Vy], beginning at OiVp^ => ORl, then every ll.ltjpg^ => 
3.0Rl. 



can be detected in these regions. The electric potential 
peak becomes wider as more particles move in from out- 
side the wake. It is these structures that will ultimately 
coalesce to fill the wake. 

2.3 Ion and Electron Phase Space 

Figures ^ and |g show the phase spaces (Y, vy) for the 
electrons and ions respectively. Note the different veloc- 
ity scales on both figures; as described earlier, the ion 
thermal velocity is less than the electron thermal veloc- 
ity {Vthe = V^Ovthi)- 

Initially, the distributions are Maxwellian with a gap 
at Y = [—1, 1]Rl- This void mimics the removal of the 
solar wind particles by the Moon. In the second phase 
space plot (t = ll.liVp^ =^ S.ORl), the electrons have 
begun to move into the wake, followed by the ions. No- 
tice that the ions have been accelerated as they move 
into the wake. This is due to the ambipolar electric field 
described earlier. The electrons continue their thermal 
expansion into the void, setting up counter-streaming 
electron beams, as seen in the fourth phase space plot 
{t = 33.3uj~^ => 9.0i?L). Now the distribution is unsta- 
ble and the two-stream instability grows. The first signs 
of the instability in ion and electron phase space can 
be seen in the fifth plot {t = UAuj-^^ ^ U.ORl). At 
Y = — 0.25-Rl there is a slight kink in the phase space 
of the ion beam travelling in the — Y direction. This 
is a direct consequence of a peak in the electric poten- 



tial generated by the two-stream instability. In the next 
phase space plot {t — bb-bcj"^ ^ 15.0Rl), more dis- 
turbances in the phase spaces of the ion beams indicate 
more disruption by the two-stream instability. Now the 
disturbances can easily be seen in the equivalent elec- 
tron phase space plot. There clearly has been a phase 
space vortex generated at Y = +0.25Rl- This is one of 
the electron phase space holes and its associated prop- 
agating ion structure, seen in the ion phase space as a 
bump on the otherwise smooth beam. In the electron 
phase space plots 6 through 10, this and other electron 
phase space holes can be seen to travel in both direc- 
tions along Y at various velocities. Their associated ion 
structures can also be tracked in the same way. One 
such electron hole has been indicated in Fig. H, having 
carefully followed it with other phase space plots. 

2.4 Two-Stream Instability in the Cen- 
tral Wake Region 

As described in previous sections, the counter-streaming 
particles set up a two-stream instability, which 
ultimately cause the wake to fill in. In this section, this 
instability is studied in more detail. 

Given a distribution of two electron beams travelling 
at ±Vo and two ion beams also travelling at ±Vb, the 
linear dispersion relation is 
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where W(,e is the electron beam plasma frequency, c^bi is 
the ion beam plasma frequency. This gives an expression 
for the frequency 
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The wavenumber of the fastest growing wave mode of the 
two-stream instability can be determined by considering 
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at which the growth rate is 
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By estimating the electron and ion densities from the 
particle density plots (Figs. g(a) and (b)) to be a third 
of the ambient solar wind density, an estimate for the 
maximum growth rate can be found 
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The growth rate can be measured by summing the 
E^ in the region of interest and taking the natural log- 
arithm. Figure || shows the result of this calculation 
over the first ^ ^Sw^g^ of the simulation. The dotted 
line shows the estimated maximum growth rate, which 
follows the increase in ln(^£''^) quite closely. The ini- 
tial increase in ln(^ E"^) is due to the ambipolar electric 
field. If the counter-streaming ion beams were absent, 
the growth rate would be adjusted to ^max ~ 0.29cL'pe 
whereas if the electrons were absent, the growth rate 
would become ^max ~ O.OBwpe- This implies that the 
phase space mixing is largely controlled by the unstable 
distribution of counter-streaming electron beams, rather 
than counter-streaming ion beams. This is contrary to 
the implications by Farrell et al. p, 02] who suggest 
the counter-streaming ion beams are the cause of the 
electrostatic instability in the wake. 

The generation of the two-stream instability produces 
some waves which travel away from the wake. The num- 
ber of the fastest growing two-stream wave mode (Eq. m 
can be substituted in the dispersion relation (Eq. O) to 
obtain the frequency and thus the velocity of this wave 
mode. 
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Figure 6: ln(X]-^Y) plotted against time (equated to 
distance behind the moon) X{Rl) for the region of the 
simulation box Y = [— 0.5, 0.5]i?L- The dotted line rep- 
resents the estimated maximum growth rate for this two- 
stream instability. 



By estimating the velocity of the counter-streaming 
electron beams, an estimate for the velocity of these 
wave can be obtained. Looking at the electron phase 
space plots, Vq w 2vthe ~ 1.8 x lO^ms"^ =4> oJts/k « 
4 X lO^ms"^. Figure ^ shows a region outside the wake 
where slight increases in electric field strength were 
found. The electric field magnitude is plotted on a log 
scale to highlight these fluctuations. The two black 
lines begin approximately where the two-stream insta- 
bility begins to develop, and travel along Y at velocities 
4.2 X lO'^ms"^ and 2.0 x lO'^ms"^. Notice the electric 
field fiuctuations begin along the fastest black line. This 
implies that the fluctuations are in some way related to 
the growth of the electron two-stream instability in the 
wake. There is more structure in these fluctuations than 
can be attributed to the two-stream instability, and this 
matter shall be revisited in Sec. III-F. 

2.5 Electron Phase Space Holes 

As described earlier (see Figs. ||(a) and (b)), double 
peaked electron enhancements about a single ion peak 
are seen to travel away from the wake at a variety of ve- 
locities. This deviation from quasineutrality is revealed 
in the charge density plot (Fig. ||(c)) and can also be 
detected in the electric field plot (Fig. §(d)). Figure g 
shows a small region of the density and electric field 
plots where two electron phase space holes with differ- 
ent velocities pass almost unaltered through each other. 
By following the trajectory of the electron holes, their 
characteristics can be measured directly from the sim- 
ulation data. The width is the distance between the 
peaks in the electric field magnitude. The density is 
the average electron density in a section of fixed length, 
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Figure 7: Electric field evolution on a log scale (Vm^^ 
shows distance along the simulation box Y{Rl) against 
time (equated to distance behind the moon) X(i?i). The 
two black lines represent a point along Y, beginning from 
where the electron two-stream instability begins to grow, 
travelling at velocities 4.2 x lO^ms"^ and 2.0 x lO^ms"^. 

about the centre of the electron hole. The amplitude is 
the height of the electric potential peak formed by the 
electron hole. Figure O shows the characteristics plotted 
against time. 

The velocity of this electron phase space hole 
increases as the simulation evolves. The method of mea- 
suring the velocity results in quite a noisy time series. 
Note the particularly noisy region at about X = — 43i?L- 
At this time, the electron hole caught up with and passed 
another electron hole. The width of this electron phase 
space hole remains almost constant (disregarding the 
noise and the interaction with the other electron hole). 
The density around the electron hole, however, clearly 
increases with time. The electron hole passes the rar- 
efaction wave at about X = —40Rl. Since the ambi- 
ent electron solar wind density is at 5cm~^, the density 
around the electron hole should level off there. This 
does happen in some of the other electron holes, but the 
one under analysis just reaches the ambient solar wind 
when the simulation is stopped. The amplitude of the 
electron phase space hole clearly decreases as the simu- 
lation evolves. At X = — 43i?i,, when the two electron 
holes pass, the amplitude here will be the sum of the 
amplitudes. The top of the peak in amplitude, is at the 
time when the electron holes coincide. 

To determine the relationships between these charac- 
teristics of the electron phase space hole, they are plot- 
ted against each other on a log- log plot (Fig. |l^). Many 
of the points are due to the interaction with the other 
electron hole, but by considering Fig. H, the approximate 
positions of these can be determined. By disregarding 
these points, definite relationships between the charac- 
teristics can be found. 

From Fig. nqa, vy oc Pg where a ~ 2.9. The other 
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Figure 8: Electron density, ion density, charge density 
and electric field magnitude showing two electron phase 
space holes passing through each other. 
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Figure 9: Characteristics of an electron phase space hole, 
plotted against time (equated to a distance behind the 
Moon) X(i?i). Shown here are (a) velocity along simu- 
lation box, Vy (lO^ms^^), (b) width, A (Ad), (c) density 
about electron hole, p^ (cm"'^), (d) amplitude, A (V). 
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Figure 10: Characteristics of an electron phase space 
hole, plotted against each other on log-log scales. Shown 
here are (a) velocity along simulation box, log^olwy) 
against density about electron hole, logio(pe), (b) ve- 
locity along simulation box, logj^pl'^^v) against ampli- 
tude, log]^Q(A), (c) density about electron hole, log]^o(pe) 
against amphtude, logio(A), (d) width, log;^o(^) against 
amplitude, log]^Q(A). 

electron phase space holes exhibit similar velocity - den- 
sity relationships, giving wyocpg. All of the other elec- 
tron phase space holes exhibit various different relation- 
ships for the other characteristics. Figure [lC|b (velocity- 
amplitude plot) displays a similar structure to the the- 
oretical and simulated electron phase space holes by 
Saeki and Gemna iQ (described later). The struc- 
ture of Fig. |lO| c shows a region of low density, occur- 
ring at the beginning of the electron phase space hole's 
evolution. The more densely populated region above 
logio(Pe) ~ 0.65 occurs as the electron hole reaches 
the ambient solar wind. With other electron holes that 
spend longer in the ambient solar wind, this stretches 
out at constant density as their amplitudes decay. 

The phase space structure of these electron phase 
space holes are shown in Fig. nTl Notice the distinct 
lack of low energy electrons in the hole. The electrons 
around the hole are trapped in the potential peak and 
are ocsillating back and forth (circling in phase space). 
The double peaked electron density structure is pro- 
duced since the electron phase space density is almost 
constant around the hole, thus resulting in the peaks 
coinciding with the hole edges and the density falling 
significantly in the centre. There is another larger elec- 
tron phase space hole, just on the edge of this plot. The 
ion phase space shows a group of ions removed from 
the bulk of the ions, coinciding with the centre of the 
electron phase space hole. This is the electron hole's 
associated ion-acoustic soliton. The beam at the top of 
the ion phase space (at vy — lO^ms"^) does not perturb 
the local fields or plasma since its contribution to the to- 
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Figure 11: Electron and ion phase space [t;Y(10®ms~^) 
against Y{Rl)] of an electron phase space hole. In order 
to see the structure, a number of phase space plots over 
a short time have been combined with the electron hole 
in the same place. 

tal charge and current density is small, however, some 
of these ions are being reflected by the positive potential 
formed by the electron hole. 

Many simulations and theoretical models have exam- 
ined electron phase space holes [HGl ISi with which the 
simulated electron holes can be compared. Saeki and 
Gemna |13] simulated an electron hole with a velocity 
less than the ion sound speed. There was a background 
of ions, which interacted with the electron hole, disrupt- 
ing it in two and generating coupled states of electron 
holes and ion-acoustic solitons. The electron holes with 
a velocity above the ion sound speed do not interact 
in the same way with the background ions. In the lu- 
nar wake simulations, the electron phase space holes are 
travelling at velocities below the ion sound speed, which 
is where the holes are coupled with the ion-acoustic soli- 
tons. Later on, the velocity has increased above the ion 
sound speed and therefore the electron holes are not dis- 
rupted further but continue moving, coupled to the ion- 
acoustic solitons, just as in Rcf. 15. Different size holes 
were considered and theoretical relationships between 
amplitude and velocity obtained. Measuring the hole 
phase space area from Fig. 11 gives a size, 5* « lOXoVthe- 
The velocity at this time is « lO^ms"^, which gives a 
Mach number, M = 4. These measurements give a the- 
oretical maximum potential of about 13V which is close 
to the measured amplitude. This implies that this model 
agrees well with the lunar wake simulations, indicating 
that these are examples of coupled states of electron 
holes and ion-acoustic solitons. Since a mass ratio of 
rrii/me = 20 has been used in the simulations, it is not 
clear whether such coupled states will be found in real- 
ity. 
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Figure 12: Electric potential about electron phase space 
hole structure (solid line) , compared to theoretical shape 
(dashed line) 

Lynov et al. |16| used a slightly modified Bernstein- 
Greene-Kruskal (BGK) scheme which assumed the fol- 
lowing form of the potential 



$(X) = <^osec\?{2X/L) 



(14) 



where the full width at half amplitude, AX = ln(l + 
\/2)L. Figure |l^ shows good agreement between the 
simulated potentials (solid line) compared to their the- 
oretical shape (dashed line) (given the amplitude and 
width) for four different times. The only disagreement 
is in the width of the potential peak. This is due to 
the measurement of the widths being the distance be- 
tween the two electric field strength peaks, rather than 
full width at half maximum. Other BGK solutions (e.g. 
$(X) = $osech4(2X/L), $(X) = $oexp(-(2X/L)2) 
and $(X) = $oexp(-(2X/L)4)) yj^y slightly different 
forms of the potential, which still closely fit the simula- 
tion results. 

2.6 Bump-on- Tail Instability Outside 
the Wake 

The electrons that pass right through the wake and 
emerge on the other side, generate a bump-on-tail dis- 
tribution which is unstable. The resulting bump-on- 
tail distribution cannot be formed from the sum of two 
Maxwellians since some of the bulk electrons have been 
moving into the void. On the beam side of the dis- 
tribution, the core Maxwellian has been truncated (see 
Fig. 0). Since the core and beam populations do not 
overlap, the appropriate dispersion relation is that ob- 
tained from the cold plasma approximation. 



^pe 



-pi 



UJ"^ [Uj-Wbf 



= 1 



(15) 



Figure 13: Dispersion relation of an electron bump- 
on-tail instability with an ion background showing real 
(sohd) and imaginary (dotted) parts (V;, w 2 x lO^ms"^, 
u]h ~ ^pel\n\ plotted over the dispersion relation pro- 
duced from a region of the simulation outside the wake 
where the waves occur. The frequencies have been nor- 
malised to the electron plasma frequency (wpe) 

where lj^ is the beam plasma frequency and V^ is the 
beam velocity. Estimates of lj^ ~ i^pe/V^ and Vf, ~ 2 x 
lO^ms"^ yields the dispersion relation shown in Fig. |l^. 
The solid line represents real values of the frequency 
and the dotted line represents the imaginary part, or 
the growth rate. 

Taking a region of the electric field where the waves 
occur, but away from any other disturbances, the dis- 
persion relation can be compared with the theoretical 
one in Fig. O. The wavelengths of the fluctuations in 
the electron density are A ~ 100m, which equates to 
a wavenumber k « 0.063m~^. These waves are purely 
electrostatic and can only be easily seen in a section 
of the simulation, related to the two-stream instability 
growth in the wake as described in Sec. III-D. This does 
not mean that the waves are only present in this section, 
just that they are amplified there. It is possible to see 
very slight fluctuations outside of this section in some 
electron density plots. 

Finally, the simulations presented in Ref. 13 were 
found to have a scaling error, outlined in Ref. 14. The 
ions had been given a charge ~ 4.5 times greater than 
they should have. This resulted in extra structures form- 
ing as described in Ref. 14. All the results presented here 
are from the corrected simulations. 

3 Conclusion 

Results from high resolution simulations of the lunar 
wake have been presented here, fully resolving the elec- 
tron phase space dynamics and leading to a more com- 
plete understanding of the ion phase space dynamics. 
Analysis of the two-stream instability generated by the 



counter-streaming particle beams from cither side of the 
wake, revealed that the instability is largely due to the 
presence of the counter-streaming electrons, rather than 
the counter- streaming ions, as previously thought. By 
examining the particle densities in detail, peaked struc- 
tures associated with the propagating electron phase 
space holes, generated by the two-stream instability, 
have been found. These structures are also evident in 
the electric field and electric potential plots. Charac- 
teristics of these electron phase space holes have been 
examined, and found to be in good comparison with pre- 
vious theoretical and computational work. An electron 
bump-on-tail instability has been found just outside the 
wake, where electrons had passed through the wake and 
created an unstable distribution on the other side. 

The variations in the IMF are expected to cause sig- 
nificant differences in the structure of the lunar wake. 
Since the IMF was at an oblique angle and varying 
in magnitude throughout the WIND flyby, comparisons 
with geometrically simple simulations such as this one 
arc difficult. Therefore, in order to simulate these con- 
ditions more closely, higher dimensional simulations are 
being undertaken. Comparisons of the simulations have 
been made and are similar to those in Ref. 13. The 
main difference is that a constant electron temperature 
was simulated. 
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